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Abstract — This paper analyzes the asymptotic exponent of 
both the weight spectrum and the stopping set size spectrum for 
a class of generalized low-density parity-check (GLDPC) codes. 
Specifically, all variable nodes (VNs) are assumed to have the 
same degree (regular VN set), while the check node (CN) set is 
assumed to be composed of a mixture of different linear block 
codes (hybrid CN set). A simple expression for the exponent 
(which is also referred to as the growth rate or the spectral shape) 
is developed. This expression is consistent with previous results, 
including the case where the normalized weight or stopping set 
size tends to zero. Furthermore, it is shown how certain symmetry 
properties of the local weight distribution at the CNs induce a 
symmetry in the overall weight spectral shape function. 



I. Introduction 

Tanner codes were introduced in [1] as a generalization of 
Gallager's low-density parity-check (LDPC) codes [2]. In the 
bipartite graph representation of a Tanner code, all variable 
nodes (VNs) have the same degree and may be interpreted 
as repetition codes with the same length (regular VN set). 
Moreover, all check nodes (CNs) are generic linear block 
codes with the same length, dimension and code book (regular 
CN set). Tanner codes with an irregular VN set (i.e., VNs 
with different degrees) or with a hybrid CN set (i.e., CNs of 
different types) are more generically referred to as generalized 
LDPC (GLDPC) codes in the literamre (e.g. [3]). Note that an 
LDPC code may be viewed as a GLDPC code where all CNs 
are single parity-check (SPC) codes. An even more general 
class of codes is represented by doubly-generaUzed LDPC 
(D-GLDPC) codes [4], where the VNs are also allowed to 
be generic linear block codes. 

In this paper, a simple formula for the asymptotic (in 
codeword length) exponent of the weight distribution of a 
GLDPC code ensemble with a regular VN set and a hybrid CN 
set is developed. As usual in the literature, this exponent will 
be referred to as the growth rate of the weight distribution or 
the weight spectral shape of the ensemble, the two expressions 
being used interchangeably throughout the paper The starting 
point for deriving the above-mentioned formula is a polyno- 
mial system solution for the spectral shape that was developed 
by the authors in [5, Theorem 1]. Here, it was shown that any 
value of the spectral shape function of an irregular D-GLDPC 
code ensemble can be calculated by solving a (4 x 4) system 
of polynomial equations, regardless of the number of different 
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VN and CN types. Additional recent works relevant to the 
subject of this paper are [6]-[10]. 

As explained in Section Ull assuming transmission over 
the binary erasure channel (BEC) and iterative decoding, the 
developed formula is also valid for the asymptotic exponent 
of the stopping set size distribution upon replacing the local 
weight enumerating function (WEF) of each CN type with an 
appropriate polynomial function. 

Symmetry properties of the growth rate of the weight 
distribution are also investigated. It is proved that the weight 
spectral shape function of a variable-regular GLDPC ensemble 
is symmetric w.rt. normalized weight a = 1/2 if the local 
WEF of each CN is a symmetric polynomial. This result 
establishes a connection between symmetry properties at a 
"microscopic" level (i.e., at the nodes of the Tanner graph) 
and symmetry of the "macroscopic" growth rate function. A 
necessary condition for symmetry of the weight spectral shape 
is also developed. 

II. Preliminary Definitions 

We consider a check-hybrid GLDPC code ensemble A^„, 
where n is the codeword length (this is equal to the number 
of VNs). All VNs are repetition codes of length q > 2, with 
input-output weight enumerating function 



S(x,2/) = l + .Ty« 



(1) 



There are Uc different CN types t E Ic = {1,2,- •• ,nc}', 
for each CN type t G Ic, we denote by ht, St and rt the 
CN dimension, length and minimum distance, respectively. We 
assume that rt > 2 for all t G Ic, and that no CN has idle bits 
(i.e., its generator matrix contains no all-zero column). The 
WEF for CN type t G Ic is given by 



u=0 



i(t) 



Here Au > denotes the number of weight-u codewords 
for CNs of type t. We denote by ut the largest u E {rt,rt + 
l,...,st} such that A^u'' > 0. 

For t e Ic, pt is the fraction of edges of the Tanner 
graph connected to type-t CNs, and the polynomial p{x) is 
defined by 

p{x) ^ - -^'-1 



tela 



If E denotes the number of edges in the Tanner graph, the 
number of CNs of type t E Ic is then given by Ept/st- 
Denoting as usual /^ p{x) dx by j p, we see that the number of 
edges in the Tanner graph is given hy E = nq and the number 



of CNs is given by m = E J p. Therefore, the fraction of CNs 
of type t e /c is given by 

Pt 



It 



stjp 



(2) 



A member of the GLDPC code ensemble Ain corresponds to 
a permutation of the E edges connecting VNs to CNs. 

In the special case where there is only one CN type (Tanner 
code ensemble), we write the WEF for this CN type as A{z) = 
J2t=o ^M-^:" = 1 + J2l.=r ^uZ^- In this case, the largest m G 
{r, r + 1, . . . , s} such that A„ > is denoted by u. 

The weight spectral shape of a GLDPC code ensemble 
sequence {A^„} is defined by 



G[a) 



lim —\o%¥.M^ 

n~>oo 77, 



\N^ 



(3) 



where E^^^ denotes the expectation operator over the ensem- 
ble A1„, iV„, denotes the number of codewords of weight w 
of a randomly chosen GLDPC code in the ensemble and the 
logarithm has base e. The limit in (|3]l assumes the inclusion 
of only those positive integers 77 for which an G Z and 
Em^ [Nan] is positive. Using standard notation, we also define 
the asymptotic relative minimum distance for A^„as a* = 
inf{a > I G{a) > 0}. The ensemble sequence is said to 
exhibit good spectral shape behavior when a* > and bad 
spectral shape behavior when a* = 0. 

Although this paper is focused on the weight spectrum, the 
results developed in Section|III]can be extended to the stopping 
set size spectrum. A stopping set of a GLDPC code may be 
defined as any subset S of the VNs such that, assuming all 
VNs in S are erased and all VNs not in S are not erased, 
every CN which is connected to S cannot recover any VN in 
iSllI A local stopping set for a CN is a subset of the local code 
bits which, if erased, is not recoverable to any extent by the 
CN. All results derived in this paper for the distance spectrum 
can be extended to the stopping set size spectrum by simply 
replacing the WEF for CN type t E Ic with its local stopping 
set enumerating function (SSEF). 

We point out that the local SSEF of a CN depends on the 
decoding algorithm used to locally recover from erasures. In 
this paper, we will consider both bounded distance (BD) and 
maximum a posteriori (MAP) CN decoding. In the former 
case, the local SSEF (BD-SSEF) is given by 



^^*\z) 



St 

E 

u=rt 



St 



(4) 



In the latter case, the local SSEF (MAP-SSEF) is given by 

$(*)(z) = l+f](/.«z" 



where 



/,(*) 



(5) 



> is the number of local stopping sets (under 



MAP decoding) of size uu 



'The concept of stopping set was first introduced in [11] in the context of 
LDPC codes. When applied to LDPC codes (i.e., all CNs are SPC codes), 
the definition of stopping set used in this paper coincides with that in [11]. 

-Denoting by Gt any generator matrix for a type-t CN, a local erasure 
pattern is a local stopping set under MAP decoding when each column of Gt 
corresponding to erased bits is linearly independent of the columns of Gt 
corresponding to the non-erased bits. 



The growth rate of the stopping set size distribution of 
the ensemble sequence {A^„} for the case of BD and MAP 
decoding at the CNs, whose definition is analogous to ([3]), will 
be denoted by G^{a) and G'$(q;), respectively. Similarly the 
asymptotic relative minimum stopping set size will be denoted 
by a'^ and a|,, respectively. 

Definition 1: Let 



^ ^ tela 



(6) 



and define the function f 



f(z) = 



E7.- 

tela 



[0, M) as 

dA<''(z 



(7) 
St for all 



AW{z) ■ 
Note that we have M = 1 if and only if ut 

tele. 

Lemma 1: The function f fulfills the following properties: 

1) It is monotonically increasing for all z > 0; 

2) f(0) = f (0) = 0; 

3) lim^^+oo f (^) = M. 

Proof: We prove the first property, as the proofs of 
the second and the third properties are straightforward. The 
derivative of f (normalized w.rt. j p) is given by 

2 



A^'\z) 



tela 



dA"'(z) d^A<"(z) 



+ z 



dz2 



dA'')(z) 
dz 



[^(*)(z)]2 



The denominator of the fraction in each term in the sum is 
strictly positive for all 2 > 0. The numerator of the fraction 
in term i G /c in the sum may be expanded as 



v=rt 
St 

z{^ uAi!h^ 



rt u=rt 



u—rt 



St 

E-^ 

u—rt 



v=rt 
St St 



4*'-"-^ + EE"("-«)4*'4*'-"-*- 



ti-i 



u^rt v^rt 



Observe that in this expression, each term in the second sum- 
mation with u = V is zero, while each (w, v) term in the second 
summation (with u > v) added to the corresponding {v,u) 
term is positive for z > 0, since u{u — v)Au A^ > z^^^~^ -\- 
v{v - u)a1*^ a1*^z"+"-i = {u~ 7;)2a1*^ a1*^z"+^'-i > and 
therefore the second summation is nonnegative for z > 0. 
Since the first summation is strictly positive for z > 0, it 
follows that V(z) > for all z > 0. D 

Note that, due to Lemma 1, the inverse of f, denoted by 
f^i : [0,M) -^ R+, is well-defined. 

III. Spectral Shape of Check-Hybrid GLDPC Codes 

We next state and prove an expression for the spectral shape 
of check-hybrid GLDPC codes. 

Theorem 1 (Spectral shape of check-hybrid GLDPC codes): 
Consider a GLDPC code ensemble with a regular VN set, 
composed of repetition codes all of length q, and a hybrid 
CN set, composed of a mixture of ric different linear block 



code types. Then, the weight spectral shape of the ensemble 
is given by 

G{a) = (1 - q)h{a) ~ qa\ogr^{a) 

+ q{fp)Y.^t\ogA'^'\r\a)) (8) 

where h{a) = — a log a— (1 — a) log(l — a) denotes the binary 
entropy function. 

Proof: In [5], a polynomial system solution for the 
spectral shape of irregular D-GLDPC codes was derived. In the 
special case of D-GLDPC codes where all VNs are repetition 
codes of length q, this is given by (special case of equation 
(8) in [5]) 

G(a) = log B{xo, yo) - a log xq 

+ g ( /p) 51 7t log A(*)(zo) + glog(l - ^) (9) 

where the values of xo, yo, zo, /3 in (|9]l are found by solving 
the (4 X 4) polynomial system 

dA('>(^0) r, 

^0 dz ^ P. 

A(*)(zo) q ' 



^E 



7t- 



(10) 



teh 



a^oyo 



1 + xoy^ 
xoyl 



q -Oi , 

/3 



and 



1 + a;o2/o 9 

zpyo _ ^ 
1 + zoyo q 



(11) 
(12) 

(13) 



Note that we are certain of the existence of a unique real 
solution to the polynomial system such that xq > 0, j/q > 0, 
zq > 0, /3 > 0, due to Hayman's formula. We solve this 
system of equations sequentially for the variables (3, zq, yo 
and Xq (respectively). First, combining (fTTT i and (fT2l i yields 



P^qa . 



(14) 



Substituting ( fT4] l into (fTOl i yields f(zo) = a which may be 
written as 

zo=f-i(a). (15) 

Using (O and ([TSll in ^^ yields 



yo = 



(l-a)f-i(a) ■ 
Finally, substituting ( fT4l i and ( fT6] l into ( fTSl i yields 

1-9 



XO = 



1 — a 



(f-^(a))' 



(16) 



(17) 



Substituting dUli, ([T5]l, (O and ([TtIi into (|9ll, and simphfying, 
leads to (Ull. D 

The expression (O holds regardless of whether the ensemble 
has good or bad spectral shape behavior. Note that, according 
to ([8]), the growth rate G{a) is well-defined only for a E 
[0, M]. This is as expected due to the following reasoning. A 
codeword of weight an naturally induces a distribution of bits 
on the Tanner graph edges, anq of which are equal to 1. Also 



note that the maximum number of ones in this distribution 
occurs when a maximum weight local codeword is activated 
for each of the 7t7Ti CNs of type t E Ic, and is thus given 
by mJ2t(£i ItUt- Hence, we have anq < mJ2tei ^t^t^ l-^-' 
a< M. " 

In Appendix U it is shown how, for small relative weight a, 
dHJ simplifies to a known expression that was derived in [12] 
for Tanner codes, and extended in [13] to irregular GLDPC 
codes. 

By considering Theorem \T\ in the special case of Tanner 
codes, we obtain the following corollary. 

Corollary 1 (Spectral shape of Tanner codes): Consider a 
Tanner code ensemble where all variable component codes 
are length-q repetition codes and where all check component 
codes are length-s codes with weight enumerating function 
A{z) = 1 + X]u=r -^uz"^- The weight spectral shape of this 
ensemble is given by 

G{a) = (1 - q)h{a) - qa \og{r\a)) + Uog A{r\a)) 

(18) 
where the function f is given by (special case of dill) 



f(z) 



zA'{z) 
sA{z) 



(19) 



and f^i : [0, M) -^ M+ is well-defined, where Af = f . 
Note that, in the special case where all CNs are SPC codes, 
( fTSl l becomes equal to the spectral shape expression for regular 
LDPC codes developed in [16, Theorem 2] for the case of 
stopping sets. Also note that, in some cases, ( fTSl l can be 
expressed analytically as f^^(a) admits an analytical form. 
An example is given in Appendix |II] 

IV. Symmetry of the Weight Spectral Shape 

Consider a GLDPC code ensemble with a regular VN set 
and a hybrid CN set. In this section, we show how a symmetry 
in the overall weight spectral shape of the ensemble is induced 
by local symmetry properties in the WEFs of the CNs. 

Definition 2: For CN type t G 4, let JJt ^ {u & N|Ai*^ > 
0}. Then, we define 

[/(**) = {u G n\ut -vEUt}. (20) 

Note that for any CN type, we always have G [/(**) and 
ut G [/(**). 

Definition 3: The WEF of CN type t E Ic is said to be 
symmetric if and only if A,g _„ = Au for all u G Ut. 
Note that if the WEF of CN type t G Ic is symmetric, then 
we have [/<**) ^ [/<*). 

Lemma 2: The WEF of CN type t G Ic is symmetric if and 
only if the all-1 codeword belongs to this code. 
The proof of Lemma |2] is omitted due to space constraints^ 

^The sufficient condition (if a linear bloclc code has the all-f codeword then 
its WEF is symmetric) is a well-loiown result in classical coding theory. On 
the other hand, we proved the necessary condition by assuming that CN type 
t G Ic has a symmetric WEF and a maximum codeword weight D < st, and 
by showing that these assumptions lead to a contradiction. As pointed out by 
one of the anonymous reviewers, the necessary condition may also be proved 
by reasoning on the expected weight of a randomly selected codeword in a 
linear block code, under the symmetry hypothesis. 



Lemma 3: The WEF of CN type t e Ic fulfills 

for all z e R+ if and only if it is symmetric. 
Proof: We have 



(21) 



AW (z-') 









where the final equality is obtained by v ~ ut — u. The 
proof is completed by observing that, if and only if A'-'^'>{z) 
is symmetric, we have [/(**) — Ut and A\^ '_^ = Ai ' for all 

veUf. ' D 

Lemma 4: The function f defined by (|7]i fulfills 



f (z) = A/ - f (z- 



(22) 



V z G M+ if and only if A^^^ (z) is symmetric for every t e /g. 
Proof: First we note that if and only if the WEF of CN 
type t G /c is symmetric, we have, differentiating ( |2TI ). 



dT" 



^ut-2 



dyl(*)(z-i) 



+ utz"'-iA(*)(2~^). 



Multiplying by z and using (l2Tl i yields 

z'-^ = -z^-^'-^^P^^U.A^Hz) 

aZ riT- i 

Then, 



dz- 



(23) 



M - f (z^i) = il/ 



E 



7t 



-l dA<''(z~^) 
dz"i 

A(*)(z-i) 






(b) 






A(*)(z) / 

u^A^^){z) + z^^-utA^^){z) 



yl(*)(z) 



f(z) 



where we have used (|6]l and ( ISTT i in (a), and ( [23T l in (b). D 
Lemma 5: The inverse function f~^ fulfills 



f-^(M-a) = — 



1 



(a) 



(24) 



V a e (0,1) if and only if A'^*^(z) is symmetric for every 

Proof: By Lemma |4] the function f fulfills ( |22] | if and 
only if A(*)(z) is symmetric for every t e 7^. By applying the 
inverse function to both sides of (|22] | and by letting f(z~^) = 
a for all z € M'*'\{0}, we obtain the statement. D 

Theorem 2 (Sufficient condition for symmetry): Consider a 
GLDPC code ensemble with a regular VN set, composed of 
repetition codes all of length q, and a hybrid CN set, composed 
of a mixture of ric different linear block codes. If A'^*)(z) is 
symmetric for each t E Ic, then the spectral shape of the 
ensemble fulfills 



G{M -a) = G{a) 



(25) 



for all ae (0,1). 



Proof: By Lemma |2] the hypothesis that A(*)(z) is 
symmetric for each t G I^ implies Uf = Sf for every t € Ic, 
and therefore M = 1. From ^ we have: 

G{M -a) = (1- q)h{M - a) - q{M - a) \ogr^{M - a) 



/p)5:7tiogA«(ri 

^ tela 



{M - a)) 



(a) 



(1 - q)h{M - a) - q{M - a) log 



f-i(a) 



+ 9 



^7aogA(') 



f-i(a) 



'=^* (1 - q)h{M -a)^qa \og{r^{a)) 

= G{a) 

where (a) follows from symmetry of A*^*^ (z) and Lemma |5] 
(b) from symmetry of A^^'^ (z), Lemma[3]and Q, and the final 
line from M = I. D 

Theorem 3 (Necessary condition for symmetry): Consider 
a GLDPC code ensemble with a regular VN set, composed 
of repetition codes all of length q, and a hybrid CN set, 
composed of a mixture of Uc different linear block codes. 
If the spectral shape fulfills G{M ~ a) = G{a) for all 
a G (0, 1), then M is a fixed point of the function 



r(a-) = 2 f 




(26) 



Proof: The proof is somewhat lengthy and we only 
sketch it due to space constraints. The first step consists of 
developing an expression for F{a) = G{M — a) — G{a). The 
obtained expression can be written as a function of z, denoted 
by J{z), by defining z = f~^(a). We must have J(z) = for 
all z e (0, +oo), and therefore J'{z) = for all z e (0, +oo). 
This latter condition must hold in particular for z = f^^(^), 
in which case we obtain the compact condition 



q logf ^ 



M 



which is equivalent to 

M = 2f 



(1 - q) log 



M 



M 



M 



2-M 



D 
Note that if A^*^ (z) is symmetric for each t G Ic then we 
have M = 1, which is always a fixed point of r(a-) defined 
in ( |26] | since it is possible to show that f (1) = ^ 



V. Examples 

Example 1 (Tanner code with (7, 4) Hamming CNs): 
Consider a rate i? = 1/7 Tanner code ensemble where all VNs 
have degree 2 and where all CNs are (7, 4) Hamming codes (it 
was shown in [14], [15] that this ensemble has good spectral 
shape behavior). The WEF of a Hamming (7, 4) CN is given 
by A{z) = 1 + 7z3 + 7z'' + z"^, while its local MAP-SSEF and 
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Fig. 1 . Spectral shapes of the Tanner code ensemble in Example [T] Solid: 
weight spectral shape (relative minimum distance: a* = 0.18650). Dashed: 
stopping set size spectral shape under MAP decoding at the CNs (relative 
minimum stopping set size: q^ = 0.11414). Dotted: stopping set size spectral 
shape under BD decoding at the CNs (relative minimum stopping set size: 
a* = 0.01025). 
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Fig. 2. Weight spectral shape of the check-hybrid GLDPC code ensemble 
in Example 12] Relative minimum distance: a* = 0.028179. 



l+7z3- 
F 212^4 



-10z^+21z^+7z*'+z' 



Iz^ + z^ respectively. 



BD-SSEF are given by $(z) = 
and *(2) = 1 + 35^3 + 352-* 
Note that we have 7\/ = - = 1 in all three cases. A plot 
of G{a), G$(a) and G^ia) obtained by implementation of 
(fTsT i is depicted in Fig. [T] We observe that A{z) satisfies the 
conditions of Theorem |2l This is reflected by the fact that 
the weight spectral shape Gijx) is symmetric with respect to 
a = 1/2. 

Example 2 (Check-hybrid ensemble): Consider a rate R = 
1/3 check-hybrid GLDPC code ensemble where all VNs 
are repetition codes of length q ^ 3 and whose CN set 
is composed of a mixture of two linear block code types 
(Ic = {1,2}). CNs of type 1 e /c are length-7 SPC codes 
with WEF A(i) (z) = [(1 + z)^ + (1 - 2)^]/2 and 71 = 0.722, 
while CNs of type 2 e h are (7, 4) codes with WEF A^^) (z) = 
l + 52^ + 7z'* + 3z^ and 72 = 0.278. The weight spectral shape 
of this ensemble, obtained from ([8]l, is depicted in Fig. |2] Note 
that for this ensemble, M = 6/7. This value is not a fixed 
point of the function F{x) defined in ( |26] | (the only fixed point 
between and 1 is .t = 0.888421). As expected, the weight 
spectral shape does not exhibit any symmetry property. 

Example 3 (Ensemble with bad spectral shape behavior): 



Fig. 3. Weight spectral shape of the Tanner code ensemble in Example [5] 



Consider a rate R = 1/5 Tanner code ensemble where 
all VNs are repetition codes of length q ^ 2 and and 
where all CNs are (5,3) linear block codes with WEF 
A{z) = 1 + 3z^ + 3z^ + z^. This ensemble is known to 
have bad spectral shape behavior (a* = 0) since we have 
A'(0)C = 6/5 > 1, where A(a;) ^ x and C = 2A2/S 
[12], [13]. A plot of the weight spectrum for this ensemble, 
obtained from ( fTSI l is depicted in Fig. [3] We observe that 
the plot of G{a) is symmetric, due to the fact that A{z) is 
symmetric (M = 1). As expected, the derivative of G{a) at 
a = is positive and hence a* = 0. 



VI. Conclusion 

A simple expression has been developed for both the weight 
and the stopping set size spectral shape of GLDPC code 
ensembles with a regular VN set and a hybrid CN set. Some 
known results (specifically, an expression for the spectral 
shape of regular LDPC codes and an asymptotic expression 
of the spectral shape for GLDPC codes as the normalized 
weight tends to zero) follow as corollaries of the developed 
formula. Moreover, symmetry properties of the spectral shape 
function have been discussed. A sufficient condition and a 
necessary condition for the weight spectral shape function to 
be symmetric have been identified. 



Appendix I 

Asymptotic Case a 







For small a, the expression dHJ reduces to a known formula 
first developed in [12] for Tanner codes, and extended in [13] 
to irregular GLDPC codes. This formula is here obtained as a 
simple corollary of Theorem [T] 

Corollary 2: In the limit where a ^- 0, the growth rate of 
the weight distribution of a GLDPC code ensemble with a 
regular VN set fulfills 



G{a) -^ [q 1 ) a log a 



^1^ log ( er fpJ2 J, AiA 



(27) 



where r denotes the smallest minimum distance of all CN 
types (i.e. r = mintg/^ r^), X^ denotes the set of CN types 



with this minimum distance (i.e. Xc = {t E Ic '■ rt ^ r}) 
and e denotes Napier's number. 

Proof: Let a = f{zo). From the definition of f given 
in d?), it is readily shown that if a ^ 0, we must have 
zo -^ 0. Next, note that in Q, each term t G Ic in the 
sum on the right-hand side is a rational polynomial in zq 
whose denominator tends to unity as a ^- 0, and whose 
numerator is dominated as a — !• by the term corresponding 
to the lowest power of zq. Therefore, as a -> 0, dTji becomes 
a -^ {J p) J2tex, ItrAi^'z^, or equivalently. 



l/r 



ZO 



•JpEtex^^tAr' 



(28) 



Using log(l + .t) ^- a; as a; ^- 0, we have that for every t E Ic 

log A(*'(zo) ~^ A^rtZo* as a ^ 0, and so 



5]7tlogA«(zo)^ ^7*4*)z5 = 



t£la 



tex. 



■Ip 



(29) 



Note that in these steps we have again used the fact that a 
polynomial expression in zq is dominated by its lowest degree 
term as zq -^ 0. Next, observe that h{a) — > —a log a as 
a ^- and therefore we obtain (using (l28b and i29[ in ([H])) 



qa 



G{a) -^{q — 1) a log a log 



•/pEtGX,7t^r 



+ 9 P 



jIp 

which coincides with ( |27] |. 
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Appendix II 
Closed Form Expressions for the Growth Rate 

It is worthwhile to note that in some cases, ( fTSI l can be 
expressed in closed form because f^^(a) can be expressed 
analytically. This is the case, for instance, for the (3, 6) regular 
LDPC ensemble, for which f (z) = a becomes ax^+bx'^+cx+ 
d = 0, where x = z^ and (a, b,c, d) = (a — 1, 15a— 10, 15a — 
5, a). This cubic equation in x may be solved by Cardano's 
method (see, e.g., [17, p. 17]; the discriminant A = p'^ + fi^ 
is negative for every a £ (0, 1), where 

Sac - b^ 9abc - 27a^d - 2b^ 

p= no ; f^ = — 



9a2 ' ^ 54a3 

The required solution is then uniquely and analytically identi- 
fied as f^^(a) = z = ^/x where x = 2-^— pcos (0/3) — j;^ > 



and 



tan 



Similarly, the weight spectral shape of a (4, 8) regular 
LDPC ensemble may be expressed in closed form through 
the solution of a quartic equation. 
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